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Abstract 

Spherically symmetric dark energy structures are investigated in the 
framework of a generalized Chaplygin gas (GCG), which has an equa- 
tion of state of the form P = —Ajp^ . We also study these in a 
modified GCG equation of state, which includes a matter term, i.e. 
P = p — A/p°. The results of the latter are then compared with 
some observational data on low-surface-brightness galaxies which are 
supposed to be dominated by dark matter. 



1. Introduction 

The generalized Chaplygin gas is one of the promising candidates to explain 
the present accelerated expansion of the universe. It is a simple model uni- 
fying dark energy and dark matter based on a perfect fluid having negative 
pressure [1]. The equation of state of the generalized Chaplygin gas is given 



where A and a are positive constants and < a < 1 and Pch and Pch are 
the generalized Chaplygin gas pressure and density, respectively. The energy 
conservation in the context of FRW cosmology yields an expression for the 
density p in terms of the scale factor a as 



The corresponding expression for the original Chaplygin gas is the special 
case a = 1 of this equation. The generalized Chaplygin gas model, like the 
original Chaplygin gas model, evolves first as dust and then as cosmological 
constant at late times. However, in its intermediate stages, it behaves as 
a mixture of a cosmological constant and a perfect fluid with an equation 
of state P — Xp. The parameter A in the generalized Chaplygin equation 
makes this model more flexible for comparison against observations. In this 
paper, we will investigate the possible existence of spherically symmetric dark 
energy structures in the framework of the generalized Chaplygin Gas (GCG). 
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2. Structure equations and solutions 

Starting with a line element 

ds^ = - dr^ - {d9^ + sin^ 9 d<P^) (3) 

and choosing 
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2GM 



(4) 



and expressing the density in terms of the enclosed mass, the conservation 
equation T^.,^ = yields 



a' 



P'+-(p + P) = (5) 



and therefore 
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Eq.(6) is a generalization of the equilibrium Euler-Poisson equation. Indeed, 
for small values of P/p, r^P, and 2GM/r, it is possible to obtain an approx- 
imation of eq. (6) as 

P' 

— + CM = . (7) 
P 

Finally, combining eq.(6) with the equation 

M' = 47rr^p , (8) 

we can get 

1 (-'7) +47rGp = . (9) 

This is a useful approximate equation to find a scaling solution for any arbi- 
trary equation of state. 



In order to determine the scaling solution for the generalized Chaplygin gas 
model density, let us assume that the density profile, we are looking for, is 

C 



p{r) 



(10) 



Here Hq is the present day Hubble constant which corresponds to the critical 
density pc, as = SnGpc/S, C and 7 being constants to be determined. 
Combining equations (9) and (10) yields 
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1 + 7(1 + a) 



(11) 



Matching the powers and the coefficients of both terms on left-hand side of 
this equation, the values of 7 and C are found to be 
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and 
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respectively. To fix A, let us consider a generalized Chaplygin gas dominated 
by the FRW model. According to this model the square of the Hubble 
parameter is given as 



dy 
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The present day critical density is 



(14) 
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choosing a = 1 at present. The constant B can also be expressed in terms 
of the matter fraction Q, at high redshift (a <^ 1) as S = {Q,pcY^°'. The 
constant A is then found to be 



A = (1 - 

4 



(16) 



and therefore, 



with 
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Based on the assumptions in order to get the approximate equation (9), this 
solution is subject to the condition 
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< 1 

(19) 

or essentially r <S -f^o"^- As can be seen from eq.(17), the spherically symmet- 
ric dark energy density can be as large as 10^ Pc, on a scale of 100 kpc. The 
metric inside the spherically symmetric Chaplygin gas is the Schwarzschild 
metric, while the background is a FRW metric. The spherically symmetric 
dark energy is embedded in the background fluid. In order to fit this spher- 
ically symmetric dark energy properly into the background, the pressure in 
both sides of the boundary must be equal [2-3], i.e. 



FRW 



P 



sew 



(20) 



which implies that the densities in the two regions must also be equal, as we 
are dealing with the same Chaplygin gas equation of state both inside and 
outside the spherically symmetric dark energy, i.e. 



Pfrw 



Pscw 



(21) 



The density profile inside the spherically symmetric system, with Schwarzschild 
metric, is M' = iirr"^ p. The outside fluid density at high redshift (a <C 1) is 
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then given as 



PFRW = ^ for a<l . (22) 



Plugging the corresponding densities into eq.(26), we obtain 

PcC _ Pc^ 



(23) 



Solving for the radius r in terms of the scale factor a and the constant B, we 
arrive at 



r{t) = B a^+- , (24) 



with the constant 



B ^H^' 



3. Discussion and Results 

In their recent paper [4], Bertolami and Paramos studied a spherically sym- 
metric dark energy structure (they called "dark energy" ) , using a polytropic 
equation of state of negative index. They argued that there are conditions 
for such objects to form, due to density fluctuations in the background of the 
generalized Chaplygin gas. According to these authors, the condition that 
helps the initial fluctuation to grow, is that the sound velocity at the surface 
of the dark star must be less than the initial expansion velocity of the dark 
star. 

According to our solution, the expansion rate is calculated to be 



r = (3 + 3q;/2) 
and the speed of sound is given by 



a3(l+a)/2 (26) 
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The ratio of the expansion rate to the speed of sound is 



— = - (2 + a) 

Cs 2^ ' ^ 2 + a) 



Vl2 + 9a > 1 . 



(29) 



Then it is true that the condition of Bertolami and Paramos [4] holds. How- 
ever, this is not sufficient to conclude that such objects exist. An additional 
point that has to be checked, is the size of the initial density contrast. The 
average density is 



m[r) 

4 ^ 
3 Trr"* 
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which means that the density contrast is 



p{a) , 



(30) 



(31) 



4 + 3q; 

Thus this density contrast is independent of time. For particular case a = 
1, the density contrast is about 0.28, which is much larger than the value 
10~^ predicted by CMB observations. Therefore, even if the expansion rate 
of the dark star is larger than the speed of sound, the density contrast is not 
sufficient to grow and form such structures. 

Generally, a non-vanishing speed of sound is a major problem for structure 
formation in all unified models of dark energy and dark matter. Such systems 
have a characteristic scale, the sonic horizon, below which the pressure frus- 
trates gravitational structure formation. Small perturbations of scales below 
this characteristic scale die off without any noticeable effect. Combining the 
continuity and the Euler-Poisson equations for the case of vanishing shear 
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and rotation, subtracting the background, and finally changing the variable 
from t to a, one gets [5] 

a 5 +-a5 - -5{l + 5) - 3 " ^ ^ (^1^ = , 

(32) 

where S — ^ _ ^ is the density contrast, Cg is the speed of sound, given by 
p 

= ^ = (33) 

^ dp p<^ ^ ^ 

for the Chaplygin gas model. 

The hnear solution for the perturbative density contrast of eq.(19), well dis- 
cussed by Fabris et al. [29], is given as 

8per{k, a) oc a"^/^ J5/14 (4 k) , (34) 

where Jv{z) is the Bessel function and k is the comoving wave number and dg 
is the sonic horizon. This ~ a for dsk -C 1, and it oscillates 

with a decaying amplitude when dsk ^ 1. This is because the non-zero speed 
of sound, and hence non-zero pressure, opposes the effect of gravity below 
that characteristic (sonic horizon) scale. This characteristic scale is, for the 
original Chaplygin gas model, given by [5] 



, _ csda _ 2 (1 - n'f a'/' ,^ , 

~ Jo am ~ 7 Q3/2 Ho ' ^ ^ 



where = ^/A/p and O = y B/{A + B) = ^/B / pcr ■ The sonic horizon 
has a radius of about 0.18 Mpc at a redshift z ~ 20. This distance is much 
larger than the size of the largest galaxy ever observed. 

It is necessary to investigate eq.(32) beyond the linear regime, in order to 
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decide whether initial perturbations can grow to the extent that they are 
powerful enough to give rise to gravitational condensation or not. Generally, 
in the nonlinear region there is a possibility, though very small, that initial 
perturbations can grow unlimited. Bilic, Lindebaum, Tupper and Viollier 
concluded in their paper [5] that unlike the linear theory, where Sr finally 
stops growing by acoustic horizon for any value of R, the perturbation ^^(a) 
grows in the nonlinear regime to infinity at a finite a, for an initial SR{adec) 
beyond a certain limit, like the case in the dust model. 

However, they also showed in the same paper that the fraction of Chaplygin 
gas which condensates, due to this infinitely large perturbation, is only 1% 
[5]. This fraction is too small, in order to conchidc that the Chaplygin gas 
favours structure formation. Being a Chaplygin gas system with a non- 
zero speed of sound (c^ = dP/dp — aA/p^ ^ 0), the spherically symmetric 
dark energy that we are discussing is governed by the general perturbation 
principle of Chaplygin gas mentioned above. This implies that structure 
formation is impossible, as the effect of gravity is opposed by the pressure. 
Thus structure formation contradicts the existence of spherically symmetric 
dark energy structures (dark stars) in the first place. 

In addition to the density contrast problem, the rotational velocity profile of 
the dark energy structures is not consistent with observation. Starting from 
the density profile determined in eq.(17), the enclosed mass can be calculated 
as 



The resulting rotational speed inside the spherically symmetrical dark energy 




(36) 
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(it could be a halo of a galaxy or even a certain spherical patch inside a 
galaxy), is then given by 



3 2 + a ^ 



-|l/2 

(Hor)'^ . (37) 



2 4 + 3a 

For the special case a — 1, the rotational velocity is proportional to r^/^. 
This relation completely contradicts the flat rotational curve observed for 
galaxies. 

However, it is possible to fix this problem by modifying the pressure and the 
density relation [6] . This can be done by combining a Chaplygin gas pressure 
with a pressure that depends on density linearly, i.e. 

P = a' p - — . (38) 

According to this modified pressure, the body would have a fiat rotational 
curve at high density, because at such densities, only the first term is dom- 
inant. The problem arises when the density is too low. In that case, the 
second term, the generalized Chaplygin gas term, is dominant and we are 
faced with the same non-fiat rotational curve problem. It is, however, possi- 
ble to determine the maximum distance or crossover point to which the first 
term dominates. This point is where the pressures are equal. P — p and 
that of the generahzed Chaplygin gas P — —A/p"' with 

PcC 

{Hor) 
4«(4 -I- 3a 



and _ 



(l - Qi+") 



2+a 

(40) 



3(2 + a)2 

are equal. The density and the rotational velocity corresponding to the 
equation of state, P — p, of the first part in the right-hand side of eq.(29) 
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are given respectively as [7] 



and Vc^(T\/2 . (41) 



The constant cr can be fixed by comparing the second equation with the value 
of the circular velocity from a real rotational curve. 

Plugging in the corresponding values of the densities and the value of the con- 
stants A and the constant cr into eq.(24) and finally solving for the crossover 
radius Tc, we obtain 



-1 



2 + a 



4{l+a) 

(42) 



3(1- 01+" 

where Vc is the rotational velocity, which varies from galaxy to galaxy. 

In order to get some numerical insight into the the comparison between the 
crossover radius and the actual radius of some LSB (low-surface-brightness) 
and HSB (high-surface-brightness) galaxies, the crossover radius is calculated 
using the above formula for a Chaplygin gas with a = 1 for both LSB and 

HSB galaxies with their observed rotational velocity. A value of 71 km/ s/Mpc 
and 0.28 are taken for the Hubble constant, Hq and f2, respectively. 

The real radius, the rotational velocity and the corresponding calculated 
value of the crossover radius for both LSB and HSB galaxies are given in 
Table 1 as de Block and McGaugh mentioned in their paper [8]. 
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Name 




-Rmax(kpc) 


-Rcross(Mpc) 


Name 




-Rmax(kpc) 


F561-1 


52 


10.1 


3.59 


DD0154 


48 


7.6 


F563-1 


111 


17.7 


11.18 


DD0168 


55 


3.4 


F563V1 


30 


7.4 


1.57 


DDO170 


66 


9.6 


F563V2 


111 


9.2 


11.18 


N55 


87 


10.2 


F564V4 


40 


-2 


2.42 


N257 


108 


9.9 


F565V2 


51 


8.4 


3.48 


N300 


97 


10.6 


F567-2 


64 


11.3 


4.90 


N801 


222 


58.7 


F568-1 


119 


14.9 


12.42 


N1003 


115 


31.3 


F568-3 


120 


16.5 


12.57 


N1530 


79 


8.3 


F568V1 


124 


19 


13.21 


N2403 


136 


19.5 


F571-8 


133 


15.6 


14.67 


N2841 


323 


81.1 


F571V1 


73 


14.6 


5.97 


N2903 


201 


24.2 


F571V2 


45 


3.7 


2.89 


N2998 


214 


46.6 


F574-1 


100 


15.4 


9.56 


N3109 


67 


8.2 


F574-2 


40 


10.7 


2.42 


N3198 


157 


29.9 


F577V1 


30 


8.9 


1.57 


N5033 


222 


35.4 


F579V1 


100 


17.3 


9.56 


N5533 


273 


74.4 


F583-1 


85 


14.6 


7.49 


N5585 


92 


9.6 


F583-4 


67 


10 


5.25 


N6503 


121 


22.2 


U0128 


131 


42.3 


14.34 


N6674 


266 


64.5 


U1230 


102 


34.7 


9.85 


N7331 


241 


36.7 


U5005 


99 


27.8 


9.42 


U2259 


90 


7.6 


U5750 


75 


21.8 


6.21 


U2885 


298 


72.5 


U5999 


155 


15.3 


18.46 









TABLE 1: The table on the left shows the rotational speed, the radius and 
the corresponding crossover radius for LSB galaxies. Similarly, the table 
on the right shows the rotational speed, the radius and the corresponding 
crossover radius for HSB galaxies [8]. 
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As can be seen from Table 1, the real radius of both the LSB and HSB 
galaxies is of order kpc. However, the corresponding crossover radius for 
each galaxy is of order Mpc. Thus this modified equation of state works as 
far as a couple of hundreds Mpc. There is a huge difference between the 
real radius and the crossover radius. This means that the point at which 
the modified equation of state fails to hold, is out of the physical extent of 
these galaxies and hence the modified Chaplygin equation of state can be 
applied, and it can adequately explain the flat rotational curves. Therefore 
this equation could be a good alternative to the Chaplygin gas equation, as 
it could possibly alleviate the rotational curve problem which the Chaplygin 
gas causes at low densities for spherically symmetric dark energy structure. 

It is clear that a must be a constant that is not characteristic of a particu- 
lar galaxy. However, the rotational velocity Vc varies from galaxy to galaxy 
and since a depends on the rotational velocity as cr = V^/-\/2, a then differs 
from galaxy to galaxy. The main problem is generalizing such models as one 
theory which governs different galaxies. 

4. Conclusions 

The unreahstic rotational velocity, the very less value of density contrast, 
together with the small probability (1%) of getting an initial density fluc- 
tuation which can result in the formation of large-scale structure exhibited 
by the spherically symmetrical dark energy structure, leads to the inevitable 
conclusion that such objects cannot exist under the context of a generalized 
Chaplygin Gas (GCG) model. 
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We then modified the GCG equation of state to P — a"^ p — A/ p°' . At 
high density the first term dominates and the resulting rotational curve is 
consistent with the flat rotational curves observed for galaxies, but at low 
density the rotational curve is unrealistic. Two types of galaxies (LSB and 
HSB) are considered to compare the real radius of these sample galaxies with 
the radius to which the modified equation of state is working for a special 
case a = \. Fortunately, the point from which this unrealistic rotational 
curve starts to dominate is completely out of the physical size of the sample 
galaxies considered. Therefore, this modified equation of state looks capable 
of explaining the spherically symmetrical dark energy rotational curve. The 
problem is as the constant a differs from galaxy to galaxy, there is no general 
equation, which holds for all galaxies. 
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